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ABSTRACT 

In this paper we investigate the combined effect of massive binary ejection from star 
clusters and a second acceleration of a massive star during a subsequent supernova 
explosion. We call this the two-step-ejection scenario. The main results are: i) Massive 
field stars produced via the two-step-ejection process can not in the vast majority 
of cases be traced back to their parent star clusters. These stars can be mistakenly 
considered as massive stars formed in isolation, ii) The expected O star fraction pro- 
duced via the two-step-ejection process is of the order of 1-4 per cent, in quantitative 
agreement with the observed fraction of candidates for isolated- 0-star formation, iii) 
Stars ejected via the two-step-ejection process can get a higher final velocity (up to 
1.5-2 times higher) than the pre-supernova velocity of the massive-star binary. 

Key words: (stars:) binaries: general, stars: formation, stars: kinematics, (stars:) 
supernovae: general 



1 INTRODUCTION 

Considering pure number counts massive stars are by far 
only a tiny minority in the stellar population of galaxies. But 
they mainly drive galactic evolution due to their dominating 
chemical and energetic feedback. Although the importance 
of massive stars for galactic astrophysics has been accepted, 
the physical circumstances of their formation are still not 
resolved, i.e. where, why, and how they form. 

It is currently strongly debated whether the formation 
of massive stars is entirely restricted to the interior of mas- 
sive star clusters or if they can form in isolation in the galac- 
tic field. 

Indeed, on the basis of a stat i stical analysis of a sample 
of galactic O stars Ide Wit et all (|2005l ') conclude that 4±2 
per cent of all O-type stars can be considered as formed 
outside a cluster environment. They further show that this 
fraction of isolated O stars is expected if the slope of the 
cluster mass function (CMF) is /? = 1.7. This assumed CMF 
slope in the low mass star cluster regime is in disagreement 
with the slope of /3 = 2 observed in the solar neighborhood 
l|Lada fc Ladall2003l l. 

The definition of an is olated O star in th e IMF- 
Monte-Carlo simulations of Ide Wit et al.l (|200 f^) is re- 
stricted to stellar e n sembl es which contain only one O star. 
iParker fc GoodwinI (|2007l ) strengthened the definition of 
an isolated O star being an O star without B-star com- 
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panions but allowed the O star to be surrounded by a 
cluster with a ma s s of <100 M©. Using this definition 
iParker fc GoodwinI (120071 ) conclude that an observed CMF 
slope o f /3 = 2 can quantitat ively explain the statistical anal- 
ysis by Ide Wit erall (|2005l) 

The analysis by Ide Wit et all (|2004 120051 ) identifies a 
small number of O-stars which are deemed to be truly iso- 
lated in the sense of not being tr aceable to an origin in a 
cluste r or OB association. However. ICvaramadze fc Bomana 
(|2008l ) reported the existence of a bow sho ck associated 
with the 0-star HD165319, which is marked in'de Wi t et al.l 
(2005) as a very likely candidate for an 0-star formed in iso- 
lation. 

In general, two main processes exist for the produc- 
tion of high- velocity O- and B-stars: i) close encounters be- 
tween binaries and single stars or binaries an d binaries can 
result in the eject i on of ma s sive st ars ( e.g. Poveda et al.l 



1 19671 : iHoffeJ 1 19831 : iMikkolal 1 1981 1 19841 : £ 

eonar dl Il99ll ) 

The ejection velocity is of the ord er of the orbital ve- 
locity of the binary (|Heggiel Il980h . Thus tighter bina- 
ries can produce larger ejection velocities, ii) The su- 
pernova explosion of one component of th e binary leads 



to a recoil of the other component (e.g. Zwickv 



1957; 



19981 : 



'Blaauw"l96l'; 'iben & Tutukov"'l997'; 'Tauris & Taken^ 
Portcgies Zwart 2000; Hoog erwerf et al. 2001). The su- 
pernova ejection scenario only requires the existence of 
massive binaries, whereas ejection rates in the dynami- 
cal ejection scenario depend on the close encounter fre- 
quency. This frequency will be increased if massive stars 
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form in compact few-body groups (jClarke 
amm- Altenburg fc KroupallioO^ . 

Various studies on the individual ejection 
processes exist. To our knowledge only one in- 

vestigation exists combining both processes 

ijGvaramadze. Gualandris fc Portegies Zwartl I2008I '). where 
the hypothesis is explored that a hyperfast pulsar can be the 
remnant of a symmetric supernova explosion of a massive 
O star dynamically ejected from a young massive star 
cluster. In this contribution we investigate for the first time 
the combination of the dynamical and supernova ejection 
process, in which a massive binary is dynamically ejected 
from a star cluster with subsequent supernova explosion 
of one binary component with recoil of the other binary 
component. We refer to this composite ejection scenario as 
the two-step-ejection process of massive stars. 

We start our investigation in Section [2] with the calcu- 
lation of the velocity spectrum of stars which are released 
with the same velocity during a supernova from binaries with 
identical ejection velocities We then derive the probability 
that stars released by a supernova from ejected binaries can 
be traced back to their parent star cluster (Section [3|, and 
discuss the maximum possible velocity which stars can get 
in the two-step-ejection process in Section [S] 



2 COMPOUND VELOCITY SPECTRUM 

Due to dynamical interactions during close encounters of 
stars, binaries can be ejected from star clusters with the ejec- 
tion velocity, Vc- If one component of the binary explodes in 
a supernova the gravitational force acting on the other com- 
ponent decreases rapidly if the expanding supernova shell 
has passed the orbit of the stellar companion releasing it 
with nearly its orbital velocity into the galactic field. De- 
pending on the configuration of the pre-supernova binary 
and the supernova details (e.g. eccentricity, supernova mass 
loss, asymmetry of the supernova) the supernova remnan t 
may still be bound to its companion. iPortegies Zwartl (|2000l ) 
calculated that between 20 and 40 per cent of the super- 
nova runaways have neutron star companions, but less than 
1 per cent are detectable as pulsars. Independent if the su- 
pernova remnant (neutron star or black hole) still remains 
bound to the other massive star or not we speak of binary 
disintegration throughout this paper. The post-supernova 
massive-star runaway is released with a release velocity, v^, 
with respect to the centre-of-mass of the pre-supernova bi- 
nary. 

The vectorial sum of both velocities, the ejection and re- 
lease velocity, is the new velocity of the released star. We call 
this velocity the compound velocity, v^. For fixed ejection ve- 
locity, Ve, and release velocity, v^, the compound velocity, Vc, 
is distributed, because the direction of the release of the star 
from the disintegrating binary during the supernova is ran- 
domly distributed in space. The corresponding spectrum of 
the compound velocity, fc{vc), defines the number of stars, 
dN{vc), which have a compound velocity, Vc, after they have 
been released from ejected binaries. 




Figure 1. Illustration of the compound velocity. A binary is 
ejected with the velocity Vc- One component is then released with 
the velocity Vr resulting in the compound velocity, Vc, given by 
the vectorial sum. 



2.1 Calculating /c(uc) 

We start the calculation of the resulting compound veloc- 
ity spectrum with the definition of the compound velocity 
through vectorial addition. 



Vc + Vr 



(1) 



The relation between the absolute values of the velocities 
follows from the cosine theorem (Fig. [1} , 



Vc='Vc~^'Vi— 2VcVi COS t 



Differentiation leads to the relation 

dVc VcVr 



de 



■ sin 6 



(2) 



(3) 



between the compound velocity and the release angle 9. 

The orientations of the binaries are randomly dis- 
tributed in space. Thus, the released stars are 47r distributed 
with respect the centre-of-mass system of the binary. If a set 
of binaries release A'^ stars isotropically the number of stars 
dN{0) released in a small angle d6 is given by the ratio of 
the area of the small circular stripe, dA{9), with the angle 
9 and the unit sphere, An (Fig. [T} , 



dN{9) _ dA{9) 

N ^ 47r ' 
where the area of the small circular stripe is 

dA{9) = 27rsin(6l)d6l. 



(4) 



(5) 



By combining these equations the number fraction of stars 
having the compound velocity Vc is 

dN{vc) _ N 



dVc 



(6) 



As the angle 6 varies from to 27r, the allowed range of the 
compound velocity can be obtained from the cosine theorem 
above and is 



V,:\ < Vc < Vc + Vr 



(7) 



The distribution function, fc{vc), of the compound ve- 
locity is normalised by 



fc{Vc) dvc 



(8) 
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Figure 2. The form of the compound velocity spectrum fc 
(eq.H. 



and is for a constant ejection velocity, Vc, and constant re- 
lease velocity, v^, 

1 



/c(Wc) 



cQ{Vc — \Vc — ■i;r|)6("o + Vr — Vc) 



2VeVr 

where the O-mapping is defined by 



e{x) 



1 ; a; > 
; a; < 



(9) 



(10) 



The form of the compound spectrum can be seen in Fig. [2l 
Compound velocities at the high speed end are preferred. 
But the distribution function, fc{vc), flattens with increas- 
ing ejection or release velocity. The distribution is symmetric 
in Vt and Vc- Note that in eq.[8]an additional factor An v^, as 
for example in the Maxwellian distribution calculated a from 
three-dimensional Gaussian distribution function, is not re- 
quired, because fc{vc) refers already to an absolute velocity 
value and is not calculated from spatial integration over a 
three-dimensional distribution function. 

If the velocities of the ejected binaries and of the re- 
leased stars are distributed according to uncorrelated distri- 
bution functions, /c(wo) and /r(«r), then the resulting dis- 
tribution of compound velocities is calculated by integration 
over both distributions, 



fciVc) 



Mvc)fv{Vr) 



VcVr 



0(Uc - \Vc - fr|) 6(Uo + - Wc) dVc dVr . (11) 



2.2 Properties of fc{vc) 

In the following we derive some properties of the compound 
velocity spectrum. 

i) The simplest case we can consider is that if one of 
the velocities, Ve or Vr, is zero. If for example Vy converges 
against zero then the velocity spectrum converges against 
the delta-distribution 



fc{Vc) = lim /i,,(Wc) = S{Vc - Va) 



(12) 



Because fc{vc) is symmetric in and Vc, the same result 
follows for Vc — >■ 0. The compound velocity is identical to 
the ejection or release velocity. 



ii) The released stars are not necessarily faster than 
the previous binaries. They can also be decelerated. The 
fraction of stars {fJ,>v^) which are accelerated, i.e. having a 
compound velocity greater than the ejection velocity is given 
by the integral 



fc{Vc) dvc 



(13) 



The evaluation of the integral leads to 

1 / 2 



4VcVr 

e{\Vc - Vr\ ~ Vc) {vr - 2VcVr)) . (14) 

The fraction of accelerated stars in dependence of the ejec- 
tion and release velocity can be seen in Fig. [S] Two cases 
can be distinguished: 

\Vc - Vr\ > Vc : fl>v„ = 1 , (15) 

and 



Vr\ <Vc : fJ->v„ 



Vr 1 

~^ 2 



(16) 



The number fraction of accelerated stars is always larger 
than 50 per cent and all stars are accelerated for > 2 Vc- 

For the case that Vc = Wr, which will be important for 
Section [5] i.e. the ejection velocity is comparable to the or- 
bital velocity of the ejected binary, the fraction of acceler- 
ated stars is 75 per cent. 

iii) We now calculate the resulting mean compound ve- 
locity, Vc, by the integral 

vc^ I ^^c/(t^c)d^;c = ^TT-r "cl,;!": , (17) 



GVcVt 



rl 



We define 

^^min = min{i;o, Wr} , Vn 

and write 

Vy — ?^min ^max 

and 



(18) 
(19) 
(20) 

and finally the mean compound velocity can be written as 

— 1 ^min '^^max 



\Vc Vy\ — '^max 



(21) 

It follows that the mean compound velocity is always greater 
than or equal to the ejection and release velocity: 



Vc > 



1 3< 



(22) 



3 ''^max 

For the case that Vc = «r, which will be important for 
section O i.e. the ejection velocity is comparable to the or- 
bital velocity of the ejected binary, the mean compound ve- 
locity is 

4 

Vc ^ - Vc (23) 



3 BACK-TRACING PROBABILITY 

If single stars or binaries are ejected from star clusters it 
is theoretically possible to calculate their orbits backward, 
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Table 1. Back-trace probability 



Vg / km s' 



Figure 3. The fraction of accelerated stars in dependence of 
the ejection and release velocity for different release velocities 
in km s~^. 
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ij; is the radius of the identification sphere around the star cluster 
from which the binary has been ejected. ? is the distance between 
the location of the disintegration of the binary and the centre of 
the identification sphere which coincides with the centre of the 
star cluster. P is the back-tracing probability in %o. 

At the location of binary disintegration (point P) the 
identification sphere appears under the angle 2a, i.e. a solid 
angle of 27r(l — cos a). Because the directions of the released 
stars are randomly distributed over the full solid angle of 
47r, the probability that the star cluster can be identified 
is given by the ratio of the solid angle of the identification 
sphere and the full solid angle, 



P = 



(24) 



The cosine can be expressed by the identification radius, 
and the disintegration distance, 4, 



= \/l - sin^ a = \ I - 



R? 



(25) 



(26) 



Then the identification probability is 

The resulting probabilities for different ratios of the identi- 
fication radius and disintegration distance are listed in Ta- 
ble [U 



Figure 4. Illustration of the back-tracing-probability. See text 
for details. 

if the force field, through which the objects have moved in 
time, is given. In such a case the star cluster, where the stars 
have their origin, can be identified. 

If the binary disintegrates due to a supernova, the re- 
leased component suffers a strong defiection of its previous 
binary orbit, and the cluster, from where it is expelled, can 
only be identified if the angle between the deflected and pre- 
vious orbit is not too large (Fig. |4]). 

A binary is ejected from a star cluster and moves along 
path 1. After the binary has moved the distance ^ from the 
centre of the cluster, O, it disintegrates at the position P and 
one component of the binary is released. The star cluster, 
from which the star has been ejected, is identified, if the 
extrapolated path of the released star intersects a sphere 
with the identification radius, _Ri, round the star cluster. 
If the released star moves along path 2, the star cluster is 
identified (extrapolated long dashed line). If the star moves 
along path 3, the parent star cluster can not been identified 
(extrapolated short dashed line) . Only stars which have new 
orbits within the gray shaded region between the dotted lines 
can be traced back to their parent star cluster. 



4 OBSERVED STATISTICS OF RUNAWAYS 
AND APPARENTLY ISOLATED O STARS 

Massive stars and massive binaries can only be ejected 
dynamically from star clusters during the early stage of 
their life. If they form withi n compact few-body c onfigu- 
rations or trapezia systems jClarke fc PringQ Il992l ). then 
the decay time scale of these few- body systems (<1 Myr, 
iPflamm-Altenburg fc Kroupall20"o^ implies only early ejec- 
tions of massive stars. If massive stars and binaries are 
formed distributed over the star cluster, then ejections can 
only occur as long as the stellar densities are high. Due to 
gas expulsion young embedded star clusters become super- 
virial a nd start to expand and the stellar density decreases 
rapidly (jKroupa. Aarseth fc Hurlevll200ll ). The time scale of 
decrease of the stellar density is comparable to the gas ex- 
pulsion time scale, of the order of ~1 Myr. Thus the time of 
flight of ejected binaries is comparable to their maximum- 
lifetime. 

Taking a lower velocity cut-off for O-sta r runa ways of 
30 km s"'^ as considered in iGies fc BoltonI i 19861 ) and a 
mean life-time of 5 Myr of O-stars implies a disintegra- 
tion distance ^ =150 pc. For a large identification radius of 
-Ri =10 pc of the star cluster, the back-tracing probability is 
1 %o (Table [!}. Lowering the cut-off velocity of the runaway- 
star definition to 10 km s~^ results in a back-tracing prob- 
ability of 1 per cent. Thus, ejected massive binaries which 
are listed in runaway O-star surveys will produce O stars 
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which can not be traced back to their parent star cluster. 
But one might expect that in this case the massive star 
can be traced back to a supernova shelL However, as sin- 
le supernova sh eUs disappear on a time scale of 0.5-1 Myr 
Chevaheij[l973 ) it might be possible to identify the parent 
supernova of the released star only in very rare cases. 

The observationally derived runaway fraction among O 
stars vari es widely in the liter ature, (see for exam ple Ta - 
ble 13 in [Cies &: Boltonlll986l or Table 5 in Stonef Il99lh . 



iGies fc BoltonI l|l986h identified 15 stars out of a sample of 
36 runaway candidates with confirmed peculiar radial veloc- 
ities >30 km s~^. Comparing with their stated total number 
of about 90 O stars within their sample space, a runaway 
fraction of 15/90= 16 per cent results. They further con- 
clude that the true runaway fraction of O stars depends 
on the adopted velocity cut-off and may lie in the range of 
10-25 per cent. i Gies fc Bohoiil (jl986) also found the binary 
fraction among runaway O stars in their sample to be about 
10 per cent. As explained above, O stars released in a su- 
pernova in an ejected massive binary result in field O stars 
which can not be traced back to their parent star clusters. 
Thus, on the basis of the O-star runaway f raction and bi- 
nary f raction of O-runaway data published bv lGies fc BoltonI 
l| 19861 ) 1-2.5 per cent of O stars can not be traced back to 
the star cluster where they have form. These O stars will 
appear to have formed in isolation, although they were born 
in an ordinary star cluster. 

The different O- star runaway fractions in the literature 
have been unified bv lStonel (|l99lh by considering true space 
frequencies: The radial velocity spectrum of O stars is de- 
composed into two different Gaussian velocity distributions, 
which correspond to two different Maxwellian space veloc- 
ity distributions. The high velocity component has a num- 
ber fraction of /h=46 per cent and a velocity distribution 
of (jh=28.2 km s~^. By transforming individual O star run- 
away studies, which are based on individual runaway defini- 
tions, to true space frequencies based on bim odality in the 
velocity distribution of O stars, IStond 11993) achieves good 
agreement between the different individual studies. 

F r om th e runaway fraction of 46 per cent derived by 
IStond (|l99ll ) it follows that a fraction of 4.6 per cent of O 
stars will apparently form in isolation, if the run away binary 
fraction of 10 per cent bv lGies fc BoltonI (| 19861 ) is used. 

Thus, the two-step-ejection process predicts a fraction 
of O stars which have f ormed apparently in isola tion in the 
range of 1-4.6 per cent. lde Wit et al.l (|2004l2005h conclude, 
based on the actually observed positions and velocities of 
O stars, that 4±2 per cent of O stars can be considered 
as candidates of massive stars formed in isolation, because 
such stars can not be traced back to a young star cluster. 
Consequently, the process of two-step ejection can quanti- 
tatively account for the proposed fraction of massive candi- 
dates formed in isolation. 



5 MAXIMUM POSSIBLE VELOCITY 

The maximum possible velocity in the two-step-ejection pro- 
cess is Umax = Vi+Vc, i.e. if the star is released in the moving 
direction of the previous binary. Binaries can be ejected from 
star clusters due to dynamical interactions. Two common 
situations are the scattering of two binaries (B-l-B) and the 



scattering of one binary and one single star (B-fS). B-l-B- 
events lead commonly to two single runaway stars and one 
tight binary. But also high velocity binaries can be produced 
in roughly 10 per cent of all cases (jMikkola 1983). 

In the B+S event high velocity binaries must be pro- 
duced due to local conservation of linear momentum. The 
ejection velocity of the single star i s typically of the order of 
the orbital velocity of the binary (Heggie"l980'), Vo- In the 
case of equal masses the ejection velocity of the binary, Vc, 
is half of the ejection velocity of the single star, Vc = ^ Vo- 
When the ejected binary disintegrates due to a supernova ex- 
plosion then the maximum possible velocity of the released 

star is Wmax = I «o. 

If on the other hand, the companion of the massive star 
is a low-mass star then the ejection velocity of the binary 
is equal to the ejection velocity of the single star, thus we 
have We = Uo . The maximum possible velocity of the released 
(less-massive companion) star after binary disintegration is 
of the order i;max = 2 It might be expected that in the 
case of binaries consisting of a high- mass and a low-mass star 
encountered by a high-mass star the low-mass star will be 
ejected and an equal-mass binary will form. To what amount 
binaries with a large mass ratio will not suffer an exchange 
has to be quantified numerically in further studies. 



6 CONCLUSION 

Various theoretical and numerical studies exist on individual 
ejection processes of massive stars, namely the dynamical 
ejection scenario and the supernova ejection scenario, but 
the combined effect of both scenarios for the distribution of 
massive stars in a galaxy has not been considered yet. In 
this paper we investigate for the first time the implications 
of the combination of the dynamical and the supernova ejec- 
tion scenario for the O-star population of the Galactic field. 
We call this combined effect the two-step-ejection process of 
massive stars. The main results are as follows: 

i) The compound velocity, Vc, which is the vectorial sum 
of the ejection velocity, v^, of the binary from the star clus- 
ter and the release velocity, v^, with which a star is released 
during a supernova, can be larger or smaller than the pre- 
vious ejection velocity. Stars can be both, accelerated and 
decelerated. 

ii) The mean compound velocity is always greater than 
each of the initial velocities, Vc and Vi (eq. [22]). 

iii) It is very unlikely that the parent star cluster of a 
massive field star produced by the two-step-ejection scenario 
can be identified. The expected number fraction of such mas- 
sive field stars which are formed apparently in isolation can 
account quantitatively for the number of candidates for iso - 
lated massive star formation derived in Ide Wit et ai]||2005h . 

iv) Massive stars which are ejected via the two-step sce- 
nario can get higher maximum space velocities (up to 1.5 
times higher for equal-mass binary components or 2 times 
higher for significantly unequal-mass companion masses) 
than can be obtained by each process (dynamical or super- 
nova ejection) individually. 
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